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An exact formula for the L2 discrepancy of the 
symmetrized Hammersley point set 

Ralph Kritzinger * 


Abstract 

The process of symmetrization is often used to construct point sets with low Lp 
discrepancy. In the current work we apply this method to the shifted Hammersley 
point set. It is known that for every shift this symmetrized point set achieves an 
Lp discrepancy of order O (\/log N/N) for p G [l,oo), which is best possible in 
the sense of results by Roth, Schmidt and Halasz. In this paper we present an 
exact formula for the L 2 discrepancy of the symmetrized Hammersley point set, 
which shows in particular that it is independent of the choice for the shift. 
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1 Introduction and statement of the result 


The local discrepancy A(q!, (d, V) of an iV-element point set V = {a^oi • • •; in the 

nnit sqnare [0,1)^ is defined as 


A(a, /3, V) = A{[0, a) x [0, /3) ,V) - Nafd 


for a, (3 & (0,1]. In this dehnition n([0, a) x [0, /3), V) is the nnmber of indices 0 < n < 
N — 1 satisfying Xn G [0,a) x [0,/3). The Lp discrepancy of a point set V in [0,1)^ is 
defined as ^ 

ipCP) = y I A(a. /3, -P) r da d/ 3 ) ' 

for p G [1, 00 ). For p —)■ cx) we obtain the notable star discrepancy. In this work we do 
not stndy this kind of discrepancy directly, bnt it should be mentioned that there is a 
remarkable asymptotic relation between the Lp discrepancy and the star discrepancy (see 
m)- The Lp discrepancy is a quantitative measure for the irregularity of distribution 
of a point set V in [0,1)^, see e.g. [H dH US]. It is also related to the worst-case 
integration error of a quasi-Monte Carlo rule, see e.g. [SIEIDE]. It is well known that 
for every p G [1, cxo) there exists a constant Cp > 0 with the following property: for the 
Lp discrepancy of any point set V consisting of N points in [0,1)^ we have 


LpiV) > Cp 


y/\ogN 

N 


( 1 ) 
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where log denotes the natural logarithm. This was hrst shown by Roth BSl for p = 2 
and hence for all p G [2, cxo) and later by Schmidt [12] for all p G (1, 2). The case p = 1 
was verified by Halasz [B]. 

Here we consider digit shifted Hammersley point sets. Let therefore m be a positive 
integer and cr = (ai, ct 2 , ..., am) G {0,1}™' a dyadic shift. We define the point set 



+ ■■■ + 


h_ Si 

2m’ 2 


■52 Sm 



where Sj = tj © aj for all j G {!,... ,m} (the operation © denotes addition modulo 
2). The point set 'Hm{<^) contains 2"^ elements. We obtain the classical Hammersley 
point set l-im with 2™ points by choosing cr = (0, 0,..., 0). Additionally, we dehne the 
m-tuple cr* = (cr^, cr^, ..., a’^) by a* = aj (Bl for all j G {1,..., m}. Then we introduce 
the symmetrized Hammersley point set as 


'^m”(<^) := 'Hm(cr) U 

This point set has elements and can be regarded as symmetrized, since 'H^™(cr) 
may also be written as the union of 'Hm(cr) with the point set 

Figure shows examples of two symmetrized Hammersley point sets. 



Figure 1: The symmetrized Hammersley point sets for z = 1, 2, where cri = 

(0,0, 0, 0, 0,0, 0, 0) and <72 = (0,1, 0,1,0,1, 0,1). The ©2 discrepancy is 0.00255571... 
in both cases. 

The concept of symmetrizing point sets plays an important role in finding point sets 
with the optimal order of Lp discrepancy in the sense of ([^. Davenport [2] used this 
method in 1956 to construct for the hrst time a two-dimensional point set with an ©2 
discrepancy of order O (^i/IoglV/A^^, and therefore showing that the lower bound ([^ 
is sharp for p = 2. For this reason, the symmetrization method we use here is often 
referred to as Davenport’s rehection principle. 

It is known that Lp{'Hm) is only of order Ci((log A^)/A^) for all p G [1, cxd) (see 
m)- However, in P, Theorem 2] it was shown with tools from harmonic analysis 


2 












(the Haar function system and the Littlewood-Paley inequality) that the symmetrized 
Hammersley point set achieves an Lp discrepancy of order O (^v^log N/n'^ for all p G 
[1, cx)) independently of the shift cr. This order is best possible in the sense of ([^. For 
the case p = 2, this result follows already from m Theorem 2] for the slightly different 
dehnition of a symmetrization of the classical Hammersley point set Urn, namely 

:= : (x, y) e Hm} ■ 

The previously mentioned results have the drawback that they do not deliver an exact 
value for the implied constant of the leading term of the L 2 discrepancy. The aim of 
this paper is to show an exact formula for the L 2 discrepancy of which gives 

not only a concrete constant, but also demonstrates that L 2 (‘H^™(cr)) solely depends 
on the number of elements N = 2™+^ and not on the shift cr whatsoever. 


Theorem 1 Let m G N and cr G {0,1}™'. Then we have 

Tn 111 1 

{2-^^L2{Wr{-))f = ^ + y + ^ - yyiiin, 
which can he displayed in terms of the number of elements N = 2™'+^ as 




1 / logiV 4 2 2 

N V 24 log 2 ^ 3 ^ iV “ ^; 


We derive the following corollary on the point set ^^“(cr) dehned as the union of 
'Hm(cr) with the point set {(x, 1 — y) : {x,y) G 'Hm(cr)} . This point set also has 2"^+^ 
elements, where some points might coincide. 

Corollary 1 Let m G N and cr G {0,1}’". Then we have with N = 2”^+^ 


L2{Wr{cT)) 


y /log^ /y\ ^ 
ivV241og2 \N) 


Proof. From |Hl Lemma 4] we have the relation 


L2{'HTi^)) - L2{Wr{<r)) 


< 


2m,+ l 


1 

N' 


Together with Theorem [T] this inequality yields the result. 


□ 


The proof of Theorem relies strongly on techniques developed and employed in 
the papers poroiDsiin!. The methods and results of [S], where the L 2 discrepancy of 
PLmitx) was computed exactly, are particularly important in order to prove the theorem. 
We comment on those results in Remark [T| Lemma and Remark The fact that we 
can write the symmetrized Hammersley point set as a union of two shifted Hammersley 
point sets allows us to employ the same techniques in this paper. The reader is invited to 
compare Theorem to the result of Kritzer and Pillichshammer as stated in Lemma 


Remark 1 Theorem [T] shows that we cannot expect a lower L 2 discrepancy by hrst 
shifting the classical Hammersley point set and then symmetrizing it. We can therefore 
simply symmetrize the classical Hammersley point set itself. This is a remarkably easy 
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construction of a point set with very low L 2 discrepancy. However, the coefficient of 
the leading term y/\ogN/N of L 2 ('H^™(cr)) is y^l/(241og 2) ^ 0.2451..., which is 
slightly higher than for the shifted Hammersley point set 'Hrn{(y) under the condition 
that the number of ones and zeros in cr is more or less balanced. In this case 
achieves an L 2 discrepancy of optimal order of magnitude in N as shown by Kritzer 
and Pillichshammer in [Hllin]- The coefficient of the leading term of L 2 {'Hrn{<^)) is then 
^5/(192 log 2) 0.1938... (see also Lemma and Remark |^. The smallest known 

leading constant is achieved for the L 2 discrepancy of digit scrambled Hammersley point 
sets in base 22 and has the value ^278629/(2811072 log 22) k, 0.1790..., as shown in 
0 - 


Remark 2 A further exact formula for the L 2 discrepancy of a symmetrized point set 
was discovered in [T]. There the authors considered symmetrized Fibonacci lattice point 
sets, whose L 2 discrepancy is also of order O {^^/\ogN/ n'^. The leading term has a 
complicated form, but numerical results indicate that the L 2 discrepancy of these point 
sets has a constant around 0.176... This would be slightly better than the result for 
digit scrambled Hammersley point sets mentioned in Remark 


2 Auxiliary results 

Throughout this paper, we call a real number a G [0,1) m-bit if it is contained in 
the set Q(2™') := {0, ..., }• Hence, a is of the form a = ^ + • • • + where 

aj e {0,1} for all j G {1,... ,m}. We also set Q*(2™') := Q(2”^)\{0}. We write Ai(q;, (]) 
for the local discrepancy of A 2 (a,/S) for the local discrepancy of and 

Asym(a,/d) for the local discrepancy of 'R^™(cr). 

The hrst lemma, which gives an exact formula for the local discrepancy of '^^(cr), 
can be derived from a result of Larcher and Pillichshammer in [TB] and was hrst stated 
explicitely in pni Lemma 1]. The second assertion in this lemma is a consequence of 
the fact that the components of all elements in 77™,(cr) are m-bit as it has already been 
pointed out in [ini Remark 3]. Here and in the following, ||a:|| := min^gz |a; — z| denotes 
the distance to the nearest integer of a real number x. 

Lemma 1 For the local discrepancy A{a,(3) ofTimif^) have 

1. A(a,/3) = ||2“/9||(-l)'^“+^(a™_„ © am+i-j{u)) for m-bit numbers a = f + 

-h and (d = ^ -\ -h (we set am+i = 0 /, where j{u) for 0 < u < m — 1 

is defined as 

fO ifu = 0, 

j{u) = < 0 if am+i-j = /3j © (Tj for j = 1 ,..., u, 

[max{j < u : am+i-j 7 ^ (dj © aj} otherwise. 

2. A(q;, 1) = 0 for m-bit a and A{a,/d) = A(a(m),/3(m)) + 2'^{a{m)l3{m) — a/d) 
for arbitrary a,/d & ( 0 , 1 ], where a{m) and /d{m) are the smallest m-bit numbers 
greater than or equal to a or /d, respectively. (For a, (d > 1 — 2“™ we choose 
a{m) = 1 and (d{m) = 1 , respectively.) 
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Lemma 2 For all a, /? G (0,1] we have Asyin(a, /9) = Ai(a, (3) + A 2 (a, (3). 

Proof. We have 

A,y^(«,/3) =A([0,«) X [0,/3),7{r(^)) 

=A([0, a) X [0, (3) , + A([0, a) x [0, {3 ), - 2"*a/? - 2"^a/3 

=Ai(a,/3) + A2{a,/3) 

for all a, (3 ^ (0,1]. □ 

Throughout the next lemma, we always write ji{u) if the function j{u) appearing in the 
hrst part of Lemmarefers to Ai{a,(3) and j 2 {u) if it refers to A2{a,(3). 

Lemma 3 Let a = ^ + ■ ■ ■ + |^ and /S = "y + ■ ■ ■ + be m-bit. 

1. For ui,U 2 G {0,..., m — 1} with Ui ^ U 2 we have 

'y (c^m—Ml © 0 'm+l—ji(ui)') (^0'm—U2 © ®m+l—^2(^2)) 2 

cgQ*(2m) 


2. For M G {0,..., m — 1} we have 

y ' iS^m—u © Q^m+ 1 —ji( m)) (o^m—n © t^m+ 1 —^2(11)) 
aeQ*( 2 ’") 

z/mG{0,1}, 

-“-1 (1 + Y.U © 1 ) 7 . + lAln © 1))) z/n G {2,... ,m - 1}. 

In the last expression, we define ■jj := (3 j 0) aj for all j G — 1}. 


Proof. We mention that Pillichshammer showed in m Lemma 2] the formula 

k 

QeQ*( 2 "‘) i=l 

for an integer 1 < A; < m — 1 and numbers Ui,... ,Uk G {0,...,m — 1} with Ui 7 ^ Uj for 
^ A i j A k, where j{u) refers to the local discrepancy of the classical Hammersley 
point set. By studying his proof, one sees that the argumentation does not change at 
all if we replace some of the j{ui) appearing in the formula by ji{ui) and the others by 
j 2 {ui), and thus we obtain the claimed identity stated in the hrst item of this lemma by 
choosing k = 2 and replacing j{ui) by jfiui) and j{u 2 ) by j 2 {u 2 )- 

We show the second item. For n = 0 we have ji{u) = 0 and j 2 {u) = 0 by dehnition 
and hence 

(c^m © Oim+l){oim © ttm+l) = J]] C^m = J]] 1 = 2 ”^ ^ = 2 ™ “ ^. 

Q-gQ*( 2 »") ai,...,«m =0 ,Om-l =0 

If M = 1, we use the fact that ji(l) and ^ 2 ( 1 ) only depend on am and write 

X] (ftm-l © am+l-ji(l)){a m— 1 © Q^m+l -i2(i)) 
aeQ*( 2 ™) 
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'y ' I y ' (c^m—1 © C^m+l—ji(1)) (oim—1 © C^m+1—j2(l)) 1 

«m=0 \ai,...,am-i=0 / 

1 

=2 (^Clm+l-ji(l)Clm+l-j2(l) © (*^«i+l-il(l) © 1) (*^™+l-i2(l) © 1)) • 

CKm—0 

We have to distinguish between the cases am = 7 i and am = 71 © 1. In the hrst case 
we obviously have ji(l) = 0 and ^ 2 ( 1 ) = 1 whereas in the second case we have ji(l) = 1 
and ^ 2 ( 1 ) = 0. We conclude 


i 

2 ^ ^ ji(l)Q^m+l—72(1) (^m+l—ji(l) © 1) (*^m+l—_ 72 ( 1 ) © 

0 

= 2 ™ ^ ^ (am © 1 ) + 2 ”^ ^ ^ (am © 1 ) 

Om= 7 l Clm =7101 

= 2 ™“ 2 ( 7 i © 1) + 2^-271 = 2™“^ = 2'"“““^ 

We turn to the case u G {2,...,m — 1}. Since ji(u) and j 2 (u) only depend on 
am+i-u, ■ ■ ■, Cim but not on , am-u, we observe that 

y ' (^m—u © Clm+1—(u)) (®m—11 © C>!.m+l—j 2 {u)') 

oeQ*( 2 ™) 

1 

y ' (o^m—u © (^m+l—ji(u}}(^m—u © Oim+l—j 2 {u)') 

ai,...,am =0 

1/1 \ 

= E E (o^m—u © ^m+1— ji (li)) (*^m—-u © Q^m+l - 32 {u)) 1 

Q^m + 1—0 yCli —-(j,=0 J 

1 

—2 

-jl{u)Cy.m+l-j2(u) + (<am+l-ji(ii) © 1) (<am+l-j2(«) © 1)) 

Q^m + 1 — li V • — 0 

U-1 1 

=2-“- E E ( ^rn-\-l—ji^rn+l—j 2 (u) © ji © 1) (*^m+l—_72(^i) © 

Jl^O Q^m+l — 

ii(«)=ii 
11-1 1 

+ 2 —E E ( C^m+l—ji{u)CX.m+l—j 2 © (o^m+l—ji{u) © l)(*am+l —^2 © ^)) 


J 2—0 — — 0 

j2{u)=j2 


=:Ti + T 2 . 


One might wonder why the sums over ji and j 2 end in m — 1 instead of u and why they 
do not coincide. The reason is that ji(u) G {0,..., m — 1} implies j 2 (u) = u and j 2 (u) G 
{0,... ,M — 1} implies ji(u) = u. This can be seen as follows: ji(u) G {0,... ,m — 1} 
implies am+i-u = 1u, because otherwise we would have ji(u) = u. But from the fact 
that ttm+i-u = 7 ii 7 ^ 7 « © 1, we immediately derive j 2 (u) = u. The other way round can 
be explained analogously. This means that the case j 2 (u) = m is actually contained in 
the sum over ji and reversely. We hnd 

©1 2 y ' ^am+lOim+l—j2{u) © (cim+1 © 1) (cim+1— 

(^m + l — u—'yu 


Oim-l=y 2 

C^m='yi 
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u-l 1 

jl=l 

(c^m+l-jittm+l -j2{u) “t" (c^m+l-ji © l)(Q^m+l-j2(n) © 1)) 

<^m+l — u—'yu 


—7^1 + 1 

=2"'-"-‘ (7„ @ 1) + x; 2J--‘ ((7J, @ 1)7„ + 7,.(7„ 0 1)) ■ 

Jl=l 


Similarly we argue that 

= 2”-“-‘7„ + 2”-“-' X; 2«-‘ ((7j, 0 1)7. + 7;,(7. 0 1)) ■ 

i2=i 


Adding Ti and T 2 completes the proof of the second item of this lemma. 
Lemma 4 Let jS be m-bit. 

1. For ui,U 2 G {0,..., m — 1} with Ui 7 "^2 we have 

om 

x: ii 2 "-))iiii 2 “>Dii = —. 

/3eQ*(2’") ^ 


For M G {0,..., m — 1} we have 


E l|2“/?f 


22m _|_ 22 w+1 

3 ■ 2”^+2 


□ 


Proof. The hrst formula follows from [TTl Lemma 3 a)] and the second one is [171 Lemma 
3 b)]. □ 

We introduce the parameter I = l{cr) := |{i G ,m} : ai = 0}|, i. e. I is the 

number of components of cr which are equal to zero. We use this notation for the rest 
of this paper. 


Lemma 5 We have 

i II Ai{a,/3)A2{a,/3) 

^ a,/3eQ*(2’") 


rnf P Im m 5 1 

9 ■ 22"^+2 ■ 


Proof. In this proof we write for the sake of simplicity A{a, /?, cr, u) := © 0 (^+ 1 -j(n), 

where we emphasize the dependence of j{u) on a, (3 and cr. With the hrst point of 
Lemma we get 

i H Ai(cr,/?)A2(cr,^) 

^ a,/3eQ*(2"*) 

1 ( 

E E l|2“‘/3||(-l)->«2l(a,,3,o-,,..) 

^ a,/3eQ*(2™) \«i=0 
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m —1 


X I i: l|2"/?||(-ir“.+M(a,/?,<T-,'t.2) 

\U2=0 


m—1 


22 r. 




Q:,/3eQ*(2"*) Vl=0 


m—1 


22 r. 


E 


a,/3eQ*(2"*) 


X E ||2“^/3||(-ir^+M(a,/3,^*,W2) 

\U 2=0 

m—1 ^ m—1 ^ 

E ||2“^/3||A(«,/?,^,«i) E ||2"^/3 ||x4(«,/3,^*,«2) 

/ 


2271 


E 


. wi =0 
X'^uj+i —0 

( , 

m—1 


a,/3eQ*(2’") 


772=0 

/ \o-u 2 + i =0 

\ / 777—1 


2271 


E 


I ui=0 
+1 —0 


m—1 


E ||2“^/?||24(a,A^,wi) E ||2“^/?||24(a,A^Xw2) 


a,/3eQ*(2’") 


1 772=0 

/ \<7'u2 + 1=1 

\ / TTl —1 


771=0 
+1 - 1 

/ 


E l|2“'/^P(a,/3,<T,Mi) E l|2“'/^P(a,/3,<T*,M2) 

72=' 

2 + 1 


1 772=0 

/ \a'uT+i =0 


\ 

) 

/ 

\ 

/ 


2271 


E 


m— 1 


\ 


I wi =0 

+ 1 — 1 


m—1 ^ 


Q,^eQ*( 2 "') 

=: — Ri + R 2 + R 3 — -R 4 . 

With the hrst part of Lemma and Lemma we obtain 


\ 772=0 


/ 


m—1 m—1 


R 2 


2271 


E E E l|2“^/5||||2 


Ui=0 U2=0 /3 gQ*(2^) 


“"'^11 E 74(a,/3,cr,Mi)74(a,/?,cr*,M2) 

QeQ*( 2 ™) 


-| rri—1 777—1 0777 1 

22777 2^ 64 ^ ’ 

Ui={j U2=0 

0 Crii2 + l~l 

In the same way we show R 3 = To calculate -Ri and i? 4 , we need to distinguish 

between the cases where Ui = U 2 and where Ui 7 ^ U 2 . This leads to 


Ri 


-I 777—1 777—1 0777 

f_ 2 ”i -2 

92777 Z^ Z^ 94 

^ 771=0 772=0 ^ 

^171 + 1—0 cr^2 + l—0 


227 


771^772 
m—1 

E E 

77=0 /3eQ*(2"i) 

^U + 1 - 0 

m—1 


77/3112 ^ A{a, (3, cr,u)A{a, (3, (T*,u) 

QgQ*(2>77) 


=s'(i-i) + T 1 : 1 : 

u =0 /3eQ*(2^) 

<Tii,+ l=0 


77/3112 ^ A{a, (3,cr,u)A{a, (3,(7*,u). 

aeQ*( 2 ™) 
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Similarly, we obtain 


1 1 

^ ^ l|2“/^ir H A{a,P,cr,u)A{a,P,cr*,u). 

^ n=0 /3eQ*(2™) aeQ*(2™) 

o'ti+i = l 

Adding i?i to R 4 yields 


— ^ Ai{a,P)A2{a,/3) 

^ a,/3eQ*(2"i) 

=- {rr? + 4/^ — 4/m — m) 

64 

1 m—1 

Y l|2“/5f II A(a,/?,cr,M)A(a,/3,<T*,M). 

^ «=0 /3eQ*(2"») aeQ*(2"*) 


Hence, our final task is to compute the last expression in the above line. We employ 
the second part of Lemma and Lemma to obtain 


m—1 


22 r. 


Y Y 

u=0 /3gQ*(2^) 


u p\\2 


H A(q!,/ 5, cr, M)A(a,/S, cr*, m) 


aeQ*(2’") 


1 


m—1 


22 r 


J2 H ||2“/3|p 

M=o /3eQ*(2’") 

2 m —1 


u—1 


+ ^ E E Il2“/5|t E 2'((7j ® 1)7« + 7j(7» ® 1)) 

m—1 2 ^ 2 m _|_ 22^+1 

'2’^+i ^ 3.2^^+^ 

-| m—1 u—1 

+ ^E2-"E2'' E l|2"/)|Pft®<r,@l)(A®<r„) 

^ 1^=2 j=l /3eQ*(2’") 


1 


m—1 71—1 


+ ^E2“"E2'' E l|2"^^IPft®<TJ)(/?„®<T„®l)-El + EJ + E3, 


u=2 j=l fiGQ*(2^) 


Finding the value of Si is a matter of straightforward calculation. We have 



For S 2 we find 

-1 m—1 u—1 1 1 

S2=lSTrE2-“E2' E E l|2"<3f. 

^*=2 j=l /3i,...,/3j_i,/9j+i,...,/3u_i=0/3„+i,...,/3m=0 

/3u=o-u(Bl 


We remark at this point that ||2“/9|p only depends on /3„+i,... ,f3m- Hence, 

1 1 o2m _j_ o2u+l 

||2“,3|p = 2 -" i: ||2”/3|p = 2-“ 1 : i|2“,3|P=2-" + 

l3u + l,...,0m=O 01,...,0m=O /3eQ*(2™) ^ 
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We arrive at 


m—1 


u—1 


2^+1 

1 

2*71+1 

m 


2““ Y 2^2““^2" 


u=2 

m—1 


Y 2 ““( 2 “ - 2 ) 2 ““^ 2 ' 


u=2 

7 


2 ^ 2 m 22 "^^+! 

3 ■ 2”^+2 

22m _j_ 22"^+! 

3 ■ 2"^+2 


96 288 9 ■ 22™+i' 

It is clear that S 3 = S 2 . Thus, after adding all the results the proof of the lemma is 
dually complete. □ 

For the proof of Theorem [T| we will also need an exact formula for the L 2 discrepancy 
of 1-irni.f^)- Such a formula was presented in P, Theorem 1], 

Lemma 6 (Kritzer and Pillichshammer) Let m eN and cr G {0,1}'”. IFe have 

’2 19m Im P I 3 m I 1 1 


m 


( 2 ’"L 2 CH„.+))) — + - + - + - + 


16 ■ 2 ^ 


■ 2 " 


+ 


4-2" 


-I- — 

I 2 ™- 


72 ■ 4”" 

, which 


Remark 3 It follows from Lemma that the optimal choice for I is 
leads to 

(2”L,(W„(o-)))* = ^+C>(l) 

(see also P Corollary 1]). This means that we achieve the optimal order of L 2 discrep¬ 
ancy for in this case. In P Theorem 1] it was shown that we achieve the optimal 

order of Lp discrepancy for all p E [1, 00 ) if and only if | 2 / — m| = 0 {y/m). 


3 Proof of Theorem [T] 


We apply Lemma to write 

(2’”+n,(Wr (<*■)))" = £ 13) f AaAI3 

= [ [ (Ai(Q;,/3))^dad/3-t- f [ {A2 {a,dadj3 
Jo Jo Jo Jo 

-h 2 y J Ai(a,/3)A2(a,/5) dad/9 

= {2^L2{n^{a))f + (2”^L2(H™(a-*)))2 

+ 2 [ [ Ai(a,/9)A2(a,/9) dad/9. (2) 

Jo Jo 

We know the values of ( 2 "^L 2 ('Hm(cr )))2 and ( 2 ”^L 2 ('Hm(cr*)))^ already from Lemma|^ 
(where in the latter case we have to insert m — l instead of I in this formula). This yields 


( 2 ™ L 2 ( H ™( a -)))2 + {2^L2inU(^*))f 


m? /2 ippi 3 X 

32^8"~”8”^W^4^ 2 "*+! 


1 

g . 22m+2 ■ 
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We split the integrals in ([^ in fonr parts: 

[ [ Ai(Q;,/9)A2(a,/3) dad/? = 

Jo Jo 


1 -2-™- rl-2- 


Ai(a, /?)A 2 (a, /?) da d/? 


f.i- 2 -m . 


+ 


+ 


,1-2-m 

f- 1 - 2 -™ 


Ai(a, /?)A 2 (a, /?) da d/? 


,1-2-m Jo 
r 


Ai(a, /?)A 2 (a, /?) da d/? 
Ai(a, /?)A 2 (a, /?) da d/? 


=:/i -|- /2 -|- /s -|- /4- 

We can calcnlate / 2 , /3 and /4 with aid of the second part of Lemma Since this 
proceeds analogonsly as in the proof of [lUl Theorem 1], we only give the resnlts. We 
have 

25 5 25 2 1 

h — h — ^ —vi: — —nr — nr—nr S" 


and 


36-2™ 9-4™ 36-4”" 3-8”" 9 ■ 16^ 

7 1 2 

— n—nr + 


6 . 4 m 9-16™ 3-8™ 

It remains to evaluate Ji. We use the second part of Lemma to obtain 


Ii = 


f-l-2-™ rl-2- 


(Ai(a(m), /?(m)) + 2™(a(m)/?(m) — a/?)) 


X (A 2 (a(m), /?(m)) + 2™'(a(m)/?(m) — a/?)) da d/? 

f-i- 2 -™ ^ 1 - 2 -™ 

/ Ai(a(m),/?(m))A 2 (a(m),/?(m)) da d/? 


" 1 - 2 -™ ri- 2 - 


+ 2 " 
+ 2 " 
+ 2 


„1_2-™ „i_2-™ 


2m 


l_2-m ._2 ... 

/ / {a(m)/3{m) — a(3)'^ da d(3 = Si + S 2 + S 3 + S 4 . 

Jo Jo 


A i(a(m), /?(m))(a(m)/?(m) — a/?) da d/3 
A 2 (a(m), /?(m))(a(m)/?(m) — a/?) da d/? 


The value of 54 can be calculated in a straightforward way and is 

1 


54 = 


;(2’"-l)^(32-2”^-25-4™-8). 


72 ■ 16”* 

The expression S 2 was computed in the proof of [91 Theorem 1] and is given by 

- 1 //(o-) m\ 


So = 2 ^ 


8 16 


Analogously, we have 


5, = 2 ™ 


_^2^+i _ I n[(T*) 


m 


8 16/’ 
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where is the number of components in cr* which are equal to zero. Since we 

obviously have /(cr*) = m — l{cr), we find S 2 + S 3 = 0. So far we have 

Ii = Si -^-(2™ - 1)2(32.2™ - 25 ■ 4”^ - 8). 

72 ■ 16^ \ ^ \ ) 

But since 


2m_i 




Ai 


a,b=l 2 ^ 



da d/3 


= II Ai(a,/3)A2(a,/3), 

^ o,/3eQ*(2™) 

we also know the value of Si from Lemma Putting all results together, we obtain the 
claimed formula in Theorem [U □ 


Acknowledgments. The author would like to thank Friedrich Pillichshammer for 
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References 

[ 1 ] D. Bilyk, V.N. Temlyakov, and R. Yu, Fibonacci sets and symmetrization in dis¬ 
crepancy theory, J. Complexity 28, No. 1 (2012) 18-36. 

[2] H. Davenport, Note on irregularities of distribution, Mathematika 3 (1956) 131- 
135. 

[3] J. Dick, F. Pillichshammer, Digital nets and sequences. Discrepancy theory and 
quasi-Monte Carlo integration, Cambridge University Press, Cambridge, 2010. 

[4] M. Drmota, R.F. Tichy, Sequences, discrepancies and applications. Lecture Notes 
in Mathematics 1651, Springer Verlag, Berlin, 1997. 

[5] H. Faure, F. Pillichshammer, G. Pirsic, and W. Ch. Schmid, L 2 -discrepancy of 
generalized two-dimensional Hammersley point sets scrambled with arbitrary per¬ 
mutations, Acta Arith. 141 (2010) 395-418. 

[6] G. Halasz, On Roth’s method in the theory of irregularities of point distributions, 
in: Recent progress in analytic number theory, Vol. 2, Academic Press, London-New 
York, 1981, pp. 79-94. 

[7] J. Horbowicz, An asymptotic relation between the extreme discrepancy and the 
Lp-discrepancy, Monatsh. Math 90, No. 4 (1980) 297-301. 

[8] A. Hinrichs, R. Kritzinger, and F. Pillichshammer, Optimal order of Lp discrepancy 
of digit shifted Hammersley point sets in dimension 2, Unif. Distrib. Theory 10 
(2015) 115-133. 

[9] P. Kritzer, F. Pillichshammer, An exact formula for the L 2 discrepancy of the 
shifted Hammersley point set, Unif. Distrib. Theory 1 (2006) 1-13. 


12 




[10] P. Kritzer, F. Pillichshammer, Point sets with low Lp discrepancy, Math. Slovaca 
57 (2007) 11-32. 

[11] L. Kuipers, H. Niederreiter, Uniform distribution of sequences, John Wiley, New 
York, 1974. 

[12] G. Larcher, F. Pillichshammer, Walsh series analysis of the L 2 discrepancy of sym- 
metrisized point sets, Monatsh. Math. 132 (2001) 1-18. 

[13] G. Larcher, F. Pillichshammer, Sums of distances to the nearest integer and the 
discrepancy of digital nets, Acta Arith. 106 (2003) 379-408. 

[14] G. Leobacher, F. Pillichshammer, Introduction to quasi-Monte Garlo integration 
and applications. Gompact Textbooks in Mathematics, Birkhauser, 2014. 

[15] J. Matousek, Geometric discrepancy. An illustrated guide. Algorithms and Gom- 
binatorics 18, Springer-Verlag, Berlin, 1999. 

[16] H. Niederreiter, Random number generation and quasi-Monte Garlo methods. Num¬ 
ber 63 in GBMS-NFS Series in Applied Mathematics, SIAM, Philadelphia, 1992. 

[17] F. Pillichshammer, On the Lp discrepancy of the Hammersley point set, Monath. 
Math. 136 (2002) 67-79. 

[18] K.F. Roth, On irregularities of distribution, Mathematika 1 (1954) 73-79. 

[19] W.M. Schmidt, Irregularities of distribution. X, in: Number Theory and Algebra, 
Academic Press, New York, 1977, pp. 311-329. 

Author’s Address: 

Ralph Kritzinger, Institut fiir Finanzmathematik und angewandte Zahlentheorie, Jo¬ 
hannes Kepler Universitat Linz, AltenbergerstraBe 69, A-4040 Linz, Austria. Email: 

r alph. kritzinger (at) j ku. at 


13 



